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Recently, a supersymmetric model of dark energy coupled to cold dark matter, the supersym- 
metron, has been proposed. In the absence of cold dark matter, the supersymmetron field converges 
to a supersymmetric minimum with a vanishing cosmological constant. When cold dark matter 
is present, the supersymmetron evolves to a matter dependent minimum where its energy density 
does not vanish and could lead to the present acceleration of the Universe. The supersymmetron 
generates a short ranged fifth force which evades gravitational tests. It could lead to observable 
signatures on structure formation due to a very strong coupling to dark matter. We investigate the 
cosmological evolution of the field, focusing on the linear perturbations and the spherical collapse 
and find that observable modifications in structure formation can indeed exist. Unfortunately, we 
find that when the growth-rate of perturbations is in agreement with observations, an additional 
cosmological constant is required to account for dark energy. In this case, effects on large scale 
structures are still present at the non-linear level which are investigated using the spherical collapse 
approach. 
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I. INTRODUCTION 

Dark energy, the component responsible for the late 
time acceleration of our Universe, is currently well de- 
scribed by a cosmological constant in the ACDM concor- 
dance model. ACDM has been very successful in explain- 
ing a large range of observations probing a vast range in 
length scales, but from a theoretical point of view the 
model suffers from the fine-tuning problem and the coin- 
cidence problem [ij. This has led to more general models 
for dark energy. Scalar field models have been particu- 
larly popular over the last decade, and are predicted to 
exist in many theories of high energy physics, like string 
theory and supergravity (see e.g. [2, 3] and references 
therein). 

However, many of the dark energy models that have 
been constructed so far suffer from problems akin to the 
ones they are trying to solve or introduce new issues 
themselves. At best they can be treated as low energy 
field theories valid well below the electron mass, corre- 
sponding to the very late phase of the Universe. Hence 
these models need to be embedded in a better defined the- 
ory whose ultra violet behaviour is under control. So far, 
no such complete scenario has been constructed. Dark 
energy models also seem to require the existence of a very 
light scalar field whose coupling to matter leads to a long 
ranged fifth force whose presence is at odds with current 
gravitational tests. Screening mechanisms (3-fToj have 
been invoked in order to alleviate this problem. Axion- 
like particles with derivative couplings to matter are also 
possible candidates (lT| . 
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On the other hand, it could well be that the dark sec- 
tor of the Universe, composed of the still undiscovered 
Cold Dark Matter (CDM) and Dark Energy (DE), could 
be described by a globally supersymmetric theory (121 - 
fl^ . In such a case the vanishingly small amount of dark 
energy which is necessary to generate the acceleration 
of the Universe could result from a small cosmological 
breaking of supersymnietry due to the non-zero CDM en- 
ergy density. Such a scenario would naturally lead to a 
close relationship between the dark energy and the CDM 
energy densities. Of course, one should also ensure that 
corrections to the globally supersymmetric scalar poten- 
tial coming from the soft supersymmetry breaking in the 
MSSM sector do not spoil the CDM-DE correspondence 
and the properties of the scalar potential in the late time 
Universe. 

Recently [isj , such a supersymmetric model of dark en- 
ergy coupled to cold dark matter, the supersymmetron, 
was proposed by two of us. In the absence of cold dark 
matter, the supersymmetron converges to a supersym- 
metric minimum with a vanishing cosmological constant. 
When cold dark matter is present, the supersymmetron 
evolves to a matter dependent minimum where its energy 
density does not vanish and can contribute to the dark 
energy budget of the Universe. 

The supersymmetron generates a short ranged fifth 
force between the CDM and the DE which evades gravi- 
tational tests, but could lead to observable signatures on 
structure formation as found in similar modified gravity 
theories [13- [11. 

In this paper we analyse the cosmological evolution 
of the supersymmetron at the background level and the 
evolution of dark matter perturbations in the linear and 
non-linear regime. The non-linear regime is studied by 
using spherical collapse. Due to the highly non-linear 
behaviour of the field during the spherical collapse we 
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are able to extract constraints on the model parame- 
ters, which are then used to constrain the background 
cosmology. The spherical collapse model has been pre- 
viously used in models with a simple Yukawa-type mod- 
ification of gravity, in the so-called f{R) / chameleon 
models [s^ |4l|, |42| , in brane- world cosmologies [i^ and 
in models which allow for dark energy fluctuations (43 |- 
[5^ . We find that a cosmological constant (CC) must be 
included in the model and that linear perturbations do 
not deviate from their ACDM counterparts. On the other 
hand, non-linear effects are significant on astrophysical 
scales. 

The outline of this paper is as follows. In Sec. ([ill '^^ 
present the supersymmetric formulation of the model, in 
Sec. (jllip we derive static solutions to the field equation 
and then in Sec. (jIVp we study the cosmological evolution 
of the supersymmetron including the cosmological back- 
ground evolution, linear perturbations and the spheri- 
cal collapse. In Sec. (|V| we revisit the original mass 
scales of the model before summarizing and concluding 
in Sec. dVll). 



given by 



II. THE SUPERSYMMETRON 

A. Supersymmetric Formulation 

In globally supersymmetric models of the scalar sec- 
tor, models are specified by their Kahler potential and 
the superpotential. With these two functions, we can 
construct the scalar potential and the kinetic term for 
the fields. For the supersymmetron we have 



if(0,0,0±,0±)- i0+r + i0- 
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where (j> is the dark energy superfield, is the CDM 
particles and A^ are some (for now) unspecified mass- 
scales. 

The kinetic term follows from 



and the scalar potential is given by the F-term 
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0n/2 
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where n — 2(/3 — a) and the mass-scales M and A are 



M 



n+4 
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Taking (p — |0|e'^, the scalar potential is seen to be min- 
imized for = —1. The angular field 6 is stabilised at 
this minimum with a mass which is always much greater 
than the gravitino mass [l5| implying that 

In the rest of this paper we will write instead of |(/)| for 
simplicity. The potential is minimized, with vanishing 
potential energy, for (f> = < 



where 

-mm = ( ^ I M 



(10) 



Due to the coupling between (/) and (j)± in the superpo- 
tential, the fermionic CDM particles acquire a scalar field 
dependent mass 



(11) 



mf{(j)) = m 1 + 



When the fermionic-CDM develops a non-vanishing num- 
ber density ricDM = (V'+V'-) in the early Universe we get 
a new contribution to the scalar potential 



VcS = Vf H PCDM, 

m 



PCDM = mncDM 



(12) 



which lifts the supersymmetric minimum and produces a 
non-zero dark energy component which can lead to the 
acceleration of our Universe. 



B. Effective 4D Model 

The effective theory for the sypersymmetron can be 
viewed as a scalar tensor theory where CDM particles 
follow geodesies of the rescaled metric g^i, = g^^A{(j)) 
where 



1 



(13) 



The effective 4D action describing the dynamics of the 
supersymmetron is given by 



9 pi 



(90)2-^^(0) 



S'cdm(-42(0)5^^; V) 



(14) 



where g is the determinant of the metric 5^1,, Mpj = 
is the reduced Planck mass and S'cdm is the dark matter 
action. If a coupling to baryons is included, the large 
mass of the supersymmetron field will ensure that this 
field would be practically invisible in local experiments. 
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C. Reparametriziation of the model parameters 



D. Supersymmetron dynamics 



In this subsection we rewrite the original mass-scales 
of the model in terms of some more intuitive physical 
quantities which will simplify our analysis. 

The coupling of the supersymmetron to dark mat- 
ter A{(j)) can be written 



^(0) + x 



where 



(15) 



(16) 



is a dimensionless parameter which parameterises the 
coupling strength of the supersymmetron to matter. 



We further introduce the density 



Poo 



(17) 



and 



Poo ^ PCDAli^ + Zoof (18) 

which is the CDM density (and the corresponding 
redshift) when the field (f> reaches the vicinity of the 
supersymmetric minimum ^min- 

The mass of the supersymmetron after having con- 
verged to the supersymmetric minimum is given by 

PaoXn 



(19) 



Constraints from supersymmetry breaking require 
TOoo 2> where 771,3/2 is the gravitino mass and is 

typically much larger than leV [la|- We will therefore 
require 

moo » C'(l)eV (20) 

which is our first constraint. The three physical param- 
eters {poo,x,moo} together with the two indices {n,/?} 
completely characterise the effective model. In the end 
we will go back and compare our results with the original 
mass-scales. 

When studying the cosmological dynamics of the 
model, it is convenient to introduce the canonically nor- 
malized field Lp via 

= n{cjy)dcjy ~> ^(0)=Ai(^-^y (21) 
In terms of Lp the potential and coupling becomes 



yf(^) = A-i 1- 



A{ip) = l + x 



Pn 



(22) 
(23) 



where (ySi, 



P{<l>n 



The field equation for ip follows from a variation of 
Eq. ([M]) with respect to p and reads 

np = Veff.v (24) 

where the effective potential is given by 

V,s{p) = Vf{p) + {A{^) - l)pcDM (25) 
The minimum pp of the effective potential is determined 

by 

Pp ) \ Pp ) Poo 

and has the approximate solution 
1 

1+1 



Pp 



PCDM ^ Po 



PCDM ^ 



(27) 



^ PCDM ^ 

A non-zero dark matter condensate is seen to lower the 
minimum from the supersymmetric minimum. The en- 
ergy density associated with the supersymmetron is 



VcsiPp) 



xpcDMjl+n) 
n 

XPCDM 



\ PCDM J 



PCDM > Po 
PCDM £ Po 



(28) 

which for small n and pcDM ^ Poo behaves like a cosmo- 
logical constant, but evolves as CDM after the field has 
converged to the supersymmetric minimum. This means 
that if the supersymmetron accounted for all dark energy, 
then acceleration would be a transient phenomenon. 

i^cs,y,y,, is given by 

\+2 



The mass of the field, = V^f 



2{n + P) (p 



2(1 



Pp 

PCDM 



Po 



Pp 



(n + 2/3) (pj^Y^ 
Pp 



(29) 
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When the field follows the minimum of the effective po- 
tential this expression simplifies to 



to; 



1 

2(n+ 



PCDM ^ Po 



(^) "'^ PCDM » Po 



(30) 



We see that the mass is always greater than the value at 
the supersymmetric minimum and from the constraint 
Eq. (|20p the mass is therefore always greater than a few 
eV's . 

The conformal coupling Eq. ([TS]) leads to a fifth-force 
(see e.g. [13, [13) between dark matter particles, which 
in the non-relativistic limit (and per unit mass) is given 
by 

rflogA((^), 



dp 



-Wp 



(31) 



This fifth-force will have an impact on structure forma- 
tion which is investigated in the following sections. 
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III. STATIC CONFIGURATIONS 

In this section we derive the static spherical symmet- 
ric solutions for the supersymmetron, which we then use 
later on when studying the spherical collapse. 

In a static spherical symmetric metric with weak grav- 
itational fields the field equation Eq. p4p reads 



(Pip 2 dip 



dr 



(32) 



We consider a spherical body of dark matter (a halo) 
with radius R and a top-hat density profile 



P = 



Pc r < R 



Pb r > R 

and impose the standard boundary conditions 

dip{r 0) dip{r — > oo) 



dr dr 

ip{r co) = ipb^ ^p[Pb) 







(33) 



(34) 
(35) 



The mass at the minimum inside (outside) the body is 
denoted by m,. (rrib). 

Outside the halo we can linearize the field equation 
around the background value ipb with the solution 



r 



r>R 



(36) 



and where the constant B is determined by matching to 
the interior solution. The interior solutions are calculated 
below for several different cases. 



A. Point particle solutions 

We first look at the point particle solution, which can 
be found by first deriving the solution for fixed R and 
then taking the limit i? — with M ~ ^R^pc fixed. 

In this limit we expect the solution inside the body to 
be a very small perturbation of the background solution 
and we can assume mbR ^1. A second order Taylor 
expansion in r gives us the solution 



A,v,Pcr^ 
6 



r<R 



Matching to the exterior solution at r = i? gives 



ip{0) = ipb ■ 



Av.PcR'' 



BR = — 
4tt 



(37) 

(38) 
(39) 



Taking the limit R ^ and using Eq. ([3T|) we find that 
the gravity plus fifth-force potential is given by 



V{r) 



GM 



(1 + 2(A ^.AfpO^e-™-'^) (40) 



which is the same as the prediction from linear per- 
turbation theorj0 as we will se later on. Contrary to 
chameleons where this type of solution holds at linear 
scales, here the large mass of the supersymmetron means 
that this solution only applies for microscopic bodies. 



B. Small overdensity 

Now we turn to the case where the size of the overden- 
sity has to be taken into account. Note that we cannot 
make the approximation rUbR 1 as the mass of the 
supersymmetron is generally very large 



111 

irib rriryo eV 



(41) 



Since we are interested in astrophysical sized overdensi- 
ties, R — O(Mpc), we will always have rribR ^ 1. 

We take ip = ip^ + 5p and Taylor expand the field 
equation inside the body around ipf) = ip{r = 0): 



d^bip 2dbip 



dr^ r dr 
which gives the solution 



ml ^ (42) 



y^^,, ( sinh(mor) 
TOor 



1 



r <R 



(43) 



Matching at r = i? and using moR, rribR 3> 1 to simplify 
the analysis, we find 



- y^b 



,ipo 



m 







sinh(moi?) cosh(mo-R) 



B : 



rrioR 



Vb - ^ t^b - y^o 



TTLbR 



(44) 
(45) 



We assume that ipo is just a small perturbation in the 
background; ipQ ~ pb ~ Sip, and expand the above ex- 
pression to first order in 6p. This leads to 



and 



Sip ■ 



B 



(Pc - Pb)A^^^ 



(Pc - P6)A' 

2ml 



(46) 



(47) 



Which gives a total force, F = '^""j^^'' on a shell close to 
the edge of the overdensity where 



Geff\ 6(A^,Mpi)^ 
G ) ^ {mbR) 



(48) 



^ By taking the Fourier transform of Eq. 1641 1 we recover Eq. 1401 1. 
see e.g. 13111 . 
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which is suppressed compared to the point-particle solu- 
tion. This solution is only valid when 



. {Pc - Pb)A^^^ 

dip ~ 2 ^ "-Pb 



Putting Pc = (1 + A)/9b and using 
that this condition reduces to 

A< 1 

i.e. a very small overdensity. 



Vb 



(49) 



mt we see 



(50) 



We have not found an explicit solution for A 1, but if 
we take A — >• in Eq. ([54| we recover Eq. (|48l) . Thus the 
two approximation agree for A ^ 1 and we will therefore 
use equation Eq. (|54p as an approximation for the fifth- 
force for all A. 

For a body of fixed size R the effective gravitational 
constant is seen to decrease with increasing pc demon- 
strating the chameleon-like behaviour and thus we re- 
cover the Newtonian regime for virialised halos, see 
Fig. ©. 



C. Large overdensity 

For large overdensities we expect screening and we 
therefore look for chameleon- like solutions [To, [25l - [29t . 
That is we assume that the field is very close to the min- 
imum, ifc = ipp{pc), of the effective potential inside the 
body and the only variations of the field are in a thin- 
shell close to the surface. Linearizing the field equation 
about ifc leads to the solution 



(p{r) ^ ipc + C 



sinh(mcr) 



r <R 



nicr 



(51) 



Due to the form of the field-equation for general /3 ^ 1 we 
cannot solve the equation in the thin-shell, but we will 
assume that this solution is valid all the way to r = i?. 
This will be the case if the shell is very thin as found in 
chameleon theories (28M30l |. and as we will see below this 
is the case for the supersymmetron when the density con- 
trast of our overdensity is large. In fact, we find that the 
supersymmetron is very similar in behaviour to strongly 
coupled chameleons as studied in (28| . 

Matching the two solutions at r = i? we obtain 



(52) 



i.e. the solution found is the critical solution where the 
field almost does not change inside the body. We can 
rewrite this equation in the standard chameleon form by 
introducing the equivalent thin-shell factor 



where $r 



AR ^ {ipb - ipc) 
R - 6/3^,$eMpV 

GAM _ {Pa-Pb)R'^ 



pi 



(53) 



R 



is the Newtonian potential 



of the overdensity. The total force on a spherical shell 
close to the surface is now F = '^"'^.^^^ where 



2/3^, (^)x(l + m,i?) (54) 



This solution is valid as long as the quadratic term in 
the Taylor expansion of Vcn.ip around ipc is suppressed 
compared to the linear term at r = i?. This condition 
turns into 



rrib 
rria 



< 1 



A » 1 



(55) 



IV. COSMOLOGICAL SUPERSYMMETRON 

In this section we discuss the cosmological evolution of 
the supersymmetron at the background level, the linear 
perturbations and the spherical collapse. 

A. Background Cosmology 

The background evolution of the supersymmetron in 
a fiat Friedmann-Lemaitre-Robertson- Walker (FLRW) 
metric 



ds^ = -df + a{tf{dx^ + dy^ + dz^) 



(56) 



is determined by the Friedman equation which in the late 
Universe reads 



Pde 



(57) 



where pb is the baryon density, pcDM the dark matter 
densiy and poE is the dark energy density. In the fol- 
lowing we will ignore the baryons and treat all matter as 
CDM. The CDM energy density is conserved implying 
that 



PCDM + 'AHpcDM — 



(58) 



The DE density is given by the sum of the energy density 
in the supersymmetron and a cosmological constant (CC) 



Pde — Pec + Pip 



(59) 



where p^ = ^ + Vf + {A{{p) — l)pcDM- We will later see 
that a non-zero CC is required to have a viable cosm olog y. 
This CC may come from supersymmetry breaking [15|. 

The field equation Eq. ([M]) in the FLRW metric 
Eq. ((56|) becomes 



(p + + Kff,^ = 



(60) 



The mass of the field is constrained by Eq. ([20)) which 
means that ^ if in the late Universe. The min- 
imum is therefore an attractor which the field fol- 
lows. Along this attractor the kinetic term is negligible 

2 



as 



(r^) ^ ^- O show the cosmo- 



logical evolution of ipp and with redshift. 
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1+Z 




l-fz 



FIG. 1. The cosmological evolution of tpp (above) and 
(below) as function of redshift. The dashed line shows the 
analytical approximation Eq. (|27p (above) and Eq. pop (be- 
low). The supersymmetron parameters are Zoo = 1.0, /3 = 1 
and n = 0.5. x and moo does not have any influence on the 
evolution of the minimum and does therefore not need to be 
specified here. 



When the supersymmetron follows the attractor we 
have 



Pv 



xpcDM 



n \f 



PCDM > P oo 



PCDM 



< 



(61) 



The equation of state along the attractor is given by 
Pv ^, Vf ^, / PCDM > Poo 



Pv 



off 







PCDM ^ poo 



(62) 



To have acceleration of the Universe without a CC we 
need to impose z^o > and n < 2. When the field 
converges to the supersymmetric minimum ~ — — >• 
and the acceleration stops. To have agreement with 
observations we need a non-zero CC, as was pointed out 
in [l^. With the inclusion of a CC, the dark energy 
equation of state is modified 



wde — 



- Pec 

Pip + Pec 



Pip > PCC 

Pip < Pec 



(63) 



_ f=1.0 
f=OJ 
f=0.1 



/ / 



E 10 

l+z 



FIG. 2. The equation of state for the supersymmetron as 
function of redshift for four different values of f = „ — : 
the fraction of dark energy in the supersymmetron to tlie total 
dark energy density today. The horizontal dotted line shows 
the analytical approximation uj^p = — The supersym- 

metron parameters are — 1.0, P = 1, n — 0.5 and x is 
fixed to give the desired / in each case. 



In Fig. ([2|) we show the dark energy equation of state 
as function of redshift and / = „ 4^0 .. " fraction 



of dark energy in the supersymmetron to the total dark 
energy density today. 

To find out how large a contribution the supersym- 
metron we can have in the energy budget of our Universe 
we will first look at the linear perturbations to get a con- 
straint on the model parameters and then apply these 
constraints to the background cosmology. 



B. Linear perturbations 

The coupling of the supersymmetron Eq. (jlSp to 
dark matter leads to a fifth-force which will influence 
the growth of the linear perturbations and structure 
formation in general. The similarity of the model with 
chameleons yields that in high density regions the 
fifth- force will be screened as shown in Sec. pilj) . 



The growth of the dark matter perturbations 6 
for sub-horizon scales are determined by 



_ SpcDM 
PCDM 

[see e.g. 



(64) 



where the effective gravitational constant is given by 



Ges{k)\ ^ , 2((l0gA)^Wpi)2 



G 



= 1 



2 1,2 



^ ^ 2(A,^Mpi)^ k 



(65) 
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where the last equahty comes from the fact that the su- 
persymmetron is very heavy compared to astrophysical 
scales and where we have assumed A{lp) — 1^1 (see 
Eq. (HHD). 

In order to have signatures on the linear perturbations 
we need the coupling strength to satisfy 2{A^tpMp\)'^ ^ 1, 
i.e. the supersymmetron must be very strongly coupled. 
It has been argued lis, 36] that an adiabatic instabil- 
ity exists in the regime, a point we will return to when 
discussing the non-linear evolution in the next section. 

At early times, pcuM ^ Poo, we find 



Gcff(fc) 
G 



1 



6 



(n+ 1)^2° 

k/a 



O.l/iMpc" 



CDM 
2 



X xlO^ 



1 + Zo 



1 



(66) 



and as the field converges to the supersymmetric mini- 
mum we obtain 



Gcff (fc) 
G 



lin 



12 



k/a 



0.1/iMpc" 



(67) 



In both cases we see that a co- moving scale of k/a = 
©(O.lft-Mpc"^) (a linear scale in GR) will experience a 
very large correction if a; ^ 1 is not satisfied. With 
X <C 1, we also have 



A{ip) -1 = X 



< X < 1 



(68) 



justifying our claim. 

To get a constraint on the model parameters we define 
^mod via 



Gcff(fcmod,^: = 0) - G 

G 



1 



(69) 



and impose fcmod > 0.1/iMpc~^ in order for the growth 
of perturbations to be in agreement with ACDM at large 
scales. With this definition we can get a constraint on 
the energy density in the supersymmetron to the total 
dark energy today. By using Eq. (p^ and Eq. we 
find 



cc 



6(1 — Ocdm) 
ing.p„10-'(l-|-2:^) 



12(l-ncDM) 



/ 0.1/t/Mpc \ 

-(^ 



l/i/Mpc 



and we must require Zoo > 10 if the supersymmetron 
is to account for a significant part of the dark energy 
budget. However, the dominating contribution to dark 
energy must be the CC as otherwise the equation of state 
Eq. (|63p reads wde ~ today and hence no acceleration. 
Thus is both cases we find that a pure CC is required to 
account for dark energy. 

Returning to the linear perturbations, we see that 
the linear effective gravitational constant is increasing 
as we go to smaller scales (large k). For k > /cmod the 
supersymmetron fifth-force is dominating over gravity 
at the linear level and to have a viable theory we 
need a large fc,„od- At non-linear scales we expect a 
chameleon-like effect to kick in and screen the fifth-force. 
We will study the non-linear effects by looking at the 
spherical collapse. 



C. Spherical collapse 

In this secion we look at the collapse of a spherical top- 
hat overdensity taking the supersymmetron fifth-force 
into account. This will allow us to extract constraints 
on the model parameters by requiring the model to agree 
with ACDM on large scales. 

The equation of motion of a spherical shell at the edge 
of the top-hat overdensity in a scalar-tensor theory with 
a fifth- force was derived in (s^. The final form of the 
equation can be understood from a simple Newtonian 
argument. In the derivation below we neglect the CC 
energy density because the Newtonian picture does not 
assign any energy density to pressure and therefore a 
Newtonian derivation cannot yield the correct contribu- 
tion (which involves contributions from pressure) from 
the CC without some ad-hoc assumptions. 

The total energy of a collapsing spherical shell of mat- 
ter is given by 



E 



1 -2 

r 



GM. 



(72) 



where V{^p) is the fifth-force potential. Neglecting shell- 
crossing so that the total energy is conserved and using 
E — Q we get Newton's law for the shell 

mshcUr = -(Fgravity(r) + F^{r)) (73) 
which can be written 



Pa 



PCBM 



m 

r{t) 



1 PCDM 



6 



1 



(70) 



-^gravity ij"^ 

1 PCBM ( Goff(r)\ 

/ sph 



6 Ml, 



G 



(74) 



If PcDM ^ Poo the energy density in the supersymmetron 
today is negligible compared to the CC. In the other 
regime Pcdm ^ P°o we find 



CC 



< 10-4(1 + Zo,)3 



(71) 



The term on the right hand side of Eq. (|74p agrees with 
the matter term found from a full derivation including 
pressure and gives the result 

rjt) _ 1 Pde _ 1 PCDM / Gcff(r)' 
r{t) ^ 3 6 I G 



(75) 



sph 
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In the following, the DE density is taken to be a pure 
CC and the effective gravitational constant is derived in 
Sec. (|III[) . see Eq. (|54p . For a small overdensity of size r 
we can write Eq. (|^5)l as 



G 



sph 



xlO' 



,36 



V r J V eV / VPCDM 



2-2/3 
n+1 



(76) 



when pcDM < Poo and 
Geff(r)' 



G 



1 



36 



/Mpc//i\ /moo A Po-- X'^JT^ 



(77) 



for pcDM 2> Poo- Note that the effective gravitational 
constant in the spherical collapse is much larger than the 
corresponding linear value. For a large-scale overdensity 
today, Mpc/ft, < r, to agree with ACDM we must require 



X < 10 



-36 



eV 



< 10 



-36 



(78) 



For such a small value of x, by looking at Eq. ()66ll67p we 
see that the linear perturbations will be indistinguishable 
from ACDM. This also means that the adiabatic insta- 
bility that might exist in these models are avoided at the 
linear level. 

By changing coordinates to y = ^ where i? = ^ we 
can write Eq. ([75)1 in the form 



y"+ (^2-^an(iV)) y' 



G 



(79) 



sph 



where a prime denotes a derivative with respect to 
N = log(a). 



The density contrast A 

ing sphere can be obtained from A = y^^ — 1 and 



iesM _ 1 of the coUaps- 

PCDM 



the mass from M 



47LR1-0 



PcDM- Early on we have 



2/ « 1 — -y as the overdensity follows the expansion, and 
by linearizing this equation we obtain the equation for 
the linear evolution of the density constrast 



where 

'Geff(r) 
G 



= 1 + 2(A <^A/pi)2(l -I- mbr)e- 



II) 



lin 



which is the same equation as for the linear perturbations 
Eq. in real space. As mentioned before 



Geff (r) - G 
G 



Geff(r)-G 
G 



< 



0{l) (82) 



sph 



on linear scales and the linear equation reduces to that of 
ACDM. This also shows that non-linear effects are very 
dominant in this model. 

The initial conditions for the numerical implementa- 
tion are taken to be the same as for ACDM: 



A, 

= 1 

3 ' 



^ A' 

3 ' ' 



A,, 



(83) 



In Fig. ([3]) we show the evolution of the radius of an over- 
density at different scales. Smaller overdensities collapse 
earlier as the fifth-force is more dominant. 

In Fig. ([4]) we plot the evolution of the effective grav- 
itational constant for the same case as Fig. ([3]). As the 
density contrast of the collapsing sphere increases, the 
chameleon mechanism kicks in and effectively shields the 
fifth-force. 

Too see more clearly the effect of the fifth-force on the 
formation of halos, we calculate the linearly-extrapolated 
density contrast for collapse today as function of the 
virial mass of the halo compared to the ACDM predic- 
tion Sc ~ 1.67, see Fig. ([S]). Low mass halos are seen 
to require a smaller linear density contrast than that of 
ACDM in order to collapse due to the fifth-force. 

With the linear collapse threshold Sc, we can predict 
the halo mass function. In the standard Press-Schechter 
approach one assumes that all regions with S > Sc 
the linear extrapolated density field collapse to form ha- 
los. The fraction of mass within halos with a given mass 
is determined by the variance of the linear density field 
smoothed over that scale. We adopt the Sheth-Tormen 
(ST) prescription [s^ for the halo mass function. The 
ST description for the co-moving number density of ha- 
los per logarithmic mass interval in the virial mass M is 
given by 



where the peak threshold 



and 



ScjM) 
a(M) 



(84) 



(85) 



A" +[2--n^{N)j A' 



lin 



(86) 



(80) We adopt the standard parameters a = 0.75 and p = 
0.3 in the following for which G = 0.322. cr(M) is the 



variance of the linear density field convolved with a top- 
hat of radius R {M = ^pi?^) 



a{Ry 



^P,ik)mkR)\'^ (87) 



where PL{k) is the linear power-spectrum, which for 
the supersymmetron is that of ACDM, and W{x) = 
(sin(a;) — x cos(x)) is the Fourier transform of the top- 
hat window function. The power-spectrum is normalized 
to (Tg = a{R = 8Mpc//i). We have chosen as = 0.8 in 
our analysis. 

The virial theorem gives us the condition for virial- 
ization of a halo. This condition reads 2T + W — 0, 
where T — j^MIi^ is the kinetic energy and W = 

— J <Pxpm{x)x ■ V^' is the potential energy. In the pres- 
ence of a fifth-force the potential 4' = $jv + A{(ppm) is 
the sum of the gravitational potential and the fifth-force 
potential. In the spherical symmetric case, we have 



W = Wn 



dA{(p) 
' dr 



dr 



where Wn = — 



3 GAr 



is the gravitational potential en- 
ergy. Performing the integration on the r.h.s. using inte- 
gration by parts we find 



R 



dA{ip) 



dr 



< 4^i?3p„(A(^(i?)) - 1) (89) 



Note that this term is usually much smaller than the 
gravitational potential energy. This can be understood 
from the chameleon thin-shell analogue: the fifth-force is 
only felt in a thin-shell close to the surface of the body 
and therefore the potential energy associated with the 
fifth-force for the whole halo is small. 

In Fig. ^ we show the ST mass function of the su- 
persymmetron relative to that of ACDM. Because the 
supersymmetron fifth-force is increasing with decreasing 
scale we recover ACDM on large scales, but see an en- 
hancement in the mass-function for low mass lialos. 



r(a) 




FIG. 3. r(a) as a function of the scale factor a for R — ^ = 
0.1, 1, 10, 100 Mpc//i (from left to right), compared with the 
behaviour for r in usual ACDM (solid). The initial density 
contrast is the same in all runs and is fixed such as to give 
collapse today for ACDM. The supersymmetron parameters 
are Zoo = 0.0, /3 = n = 1, a; = 10"''^ and moo = lO^eV. 
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FIG. 4. The strength of the fifth-force at the surface of the 
spherical overdensity during the collapse as a function of the 
scale factor a for R = ^ = 0.1, 1, 10, 100 Mpc//i (from top to 
bottom). The parameters are the same as in Fig. (1). 



V. MASS SCALES 

Having found a range for our model parameters which 
gives predictions that are in agreement with current 
observations, we will now analyse how these constraints 
affect the original mass-scales Ai of the model. 



We start by defining x' = lO^^x, m' — 
z' = {1 + z^)^. We can now rewrite Eq. ([IB] 



ffToV 

as 



Ho 



(90) 



and Eq. (H?!) as 



From Eq. pO|) we find 

U^pJ 
By using Eq. ([7]) we get 



10" 



-160-55n^m+l^/n^/ 



and lii 



from which we find 



(92) 



AoJ \MpJ 



(94) 
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FIG. 5. The ACDM-lineary-extrapolated critical density con- 
trast for collapse for the supersymmetron as a function of 
the halo mass. The dashed line shows the ACDM prediction 
Sc « 1.67. The supersymmetron parameters are Zoo = 0.0, 
13 ^n=l,x = 10""=' and = 10^ eV. 



U 
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c 




FIG. 6. The fractional difference in the supersymmetron mass 
function compared to ACDM at z = 0. The symmetron 
parameters are rrioo ~ lO^eV, Zoo = 0.0, x = 10"^'^ and 
P = n = 1. The supersymmetron converges to ACDM for 
large halo masses. 



The simplest case to analyze is /? = 1 for which the 
scale Ai vanishes from the theory. We find 



A. 



(96) 



i.e. A2 needs to be between the current Hubble scale 
and the dark energy scale. This scale can be elevated by 
increasing the redshift z^o ^ 1, but we typically need 
a redshift in the very early Universe to reach super- TeV 
scales. For Aq we find 



An 



10 



20 



(a;'™')" 



(97) 



In the general case /3 7^ 1 we see from Eq. (|93p that taking 
/3 < 1 together with Ai ^ Ao,A2 can serve to increase 



the other two scales. For example /3 
gives 



i and Ai sa Mpi 



A2 « Ao « lO^^ix'z'pHo 



(98) 



which is around the dark energy scale. There seems to be 
no unflne-tuned way of bringing these mass-scales up to 
typical particle physics scales if we want the cosmological 
symmetry breaking to be close to the present era. For 
example to have A2 ~' TeV when (3 = 1 and w 1 then 
Eq. shows that we need z' w 10^™ which translates 
to Poo ~ (lOi^GeV)^ 

Finally and from a field theoretic point of view, the 
scale Ai has a different status from Ao,2- The former ap- 
pears in the Kahler potential as a suppression of scales 
for higher dimensional operators and signals the typical 
scales above which the effective field theory description 
breaks down. On the contrary, Ao,2 appear in the su- 
perpotential and are protected by non-renormalisation 
theorems. Hence we expect that Ai should be sensitive 
to high energy physics and represents the effective cut-off 
of the theory. On the other hand, Ao,2 may be already 
present at very high energy even if these scales are very 
low. Of course, this does not provide an explanation for 
the discrepancy of scales between Ai and A0.2 which is 
not natural. 



VI. DISCUSSION AND CONCLUSION 

We have studied the cosmological evolution of the su- 
persymmetron and its possible effects on structure for- 
mation. Requiring that linear perturbations are in agree- 
ment with ACDM on large scales, we find that the energy 
density in the supersymmetron is negligible compared to 
the dark matter density and a pure cosmological constant 
must be introduced to play the role of dark energy. 

The non-linear evolution of the model was also inves- 
tigated by using the spherical collapse model. Spheri- 
cally symmetric solutions to the field equation have been 
derived and used to predict the fifth-force effects on a 
collapsing halo. The effective gravitational constant at 
the edge of a spherical overdensity has been found to be 
much larger than the linear prediction due to the highly 
non-linear properties of the model. The model parameter 
must be tuned such that the spherical collapse is under 
control on large scales. This implies that linear perturba- 
tions reduce to that of ACDM. On non-linear scales the 
model then predicts a faster collapse than that of ACDM. 
In particular we find that the supersymmetron predicts 
an excess of small mass halos compared to ACDM. How- 
ever for this to be the case, one or more of the mass-scales 
in the theory must be fine-tuned. 

On very small scales, i.e. galaxy-scales, the matter 
density is large enough to effectively screen the fifth-force 
via the chameleon mechanism. This non-linear regime 
could in principle be probed using N-body simulations. 
However, due to the enormous mass of the field this poses 
a severe challenge for existing methods. 
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